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We analyze the analog of the Rayleigh instability in the field of gravity for the superfluid-crystal 
''He interface provided that the heavier *He crystal phase occupies the half-space over the lighter 
superfluid phase. The conditions and the onset of the gravitational instability are different in kind 
above and below the roughening transition temperature when the crystal *He surface is in the rough 
or in the smooth faceted state, respectively. In the rough state of the surface the gravitational 
instability is similar to the classical case of the fluid-fluid interface. In contrast, in the case of 
the crystal faceted surface the onset of the gravitational instability is associated with surmounting 
some potential barrier. The potential barrier results from nonzero magnitude of the linear facet 
step energy. The size and the tilting angle of the crystal facet are also important parameters for 
developing the instability. The initial stage of the instability can be described as a generation of 
crystallization waves at the superfluid-crystal interface. The experiments which may concern the 
gravitational instability of the superfluid-crystal *IIe interface are discussed. 

PACS numbers: 67.80.bf, 47.20.Ma 



I. INTRODUCTION 

The Rayleigh instability is a well-known instability of 
the interface between two liquids in the field of gravity 
when the heavier liquid is placed above the lighter one. 
The similar effect is difficult to observe in solids since 
any displacement of a piece of the solid body is impeded 
due to appearing elastic stresses. However, there exists 
an elastic stress-free possibility of changing the shape 
of a crystal, namely, remelting in the hydrostatic pres- 
sure gradient. A "^He crystal in contact with its super- 
fluid phase could be one of most promising objects to 
observe the manifestation of the Rayleigh instability for 
the liquid-solid interface. In fact, solid ''He density is 
larger than that of the liquid phase. The growth rate 
of atomically rough crystal surfaces increases drastically 
with the lowering of the temperature, and the change of 
the crystal shape can occur in very short time intervals of 
about 1 s. Unfortunately, in most of crystal growth ex- 
periments Tj a ^He crystal appears either at the bottom 
of an experimental cell or drops on the same bottom later 
if a crystal seed nucleates first at the wall. Thus, prepar- 
ing the necessary configuration with the heavier crystal 
above its lighter liquid phase in order to have the initial 
condition for the development of the classical Rayleigh 
instability comes across a difficulty. 

Nevertheless, one experiment, in which the necessary 
liquid-crystal configuration against the direction of grav- 
ity is arranged, has been reported 0. A single hep ''He 
crystal is grown at the bottom of an experimental cell at 
the temperature of about 1.1 K lying between the first 
and second roughening transitions. After the grown ^He 
crystal occupies the lower half of the optical cell, the 
cell is rotated mechanically through 180° so that the ^He 
crystal proves to be above the superfluid liquid. Then 
the crystal phase starts to melt, descending along the 
cell walls. In its turn, a single finger of the superfluid 
phase moves in the upward direction at the center of the 
cell. Eventually, the crystal again occupies the lower half 




FIG. 1: The "He crystal grown at 0.92 K between two rough- 
ening transitions as in experiment [31 ■ The diameter of the 
visible margin is 12 mm. 



of the cell. Unfortunately, the development of the insta- 
bility observed is described qualitatively and the results 
of visual observations are illustrated by the schematic 
figure alone. The main conclusion of the work is that 
the spatial scale associated with the development of the 
instability is about 1 cm. 

This result [2] seems to conflict with the stability ob- 
served for the lower surface of the '*He crystals grown 
on the needle point at the center of a cell (Fig. [T]). As 
is seen from the figure, the lower crystal facet is under 
conditions appropriate for the development of the grav- 
itational Rayleigh instability. The single distinction is 
likely to be associated with the state of the crystal ^He 
surface. Whereas in experiments [3| one deals with the 
rough state of the crystal ^He surface, in our experiments 
!3] one observes the smooth faceted state of the surface in 
the form of a hep crystal basal facet above the superfluid 
phase. 

In this paper we consider an instability of the crystal 
''He surface above the superfluid phase in the field of 
gravity. We analyze and compare two possible states of 
the crystal ^He surface such as the rough and smooth 
faceted ones. 
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II. INSTABILITY OF THE CRYSTAL SURFACE 
IN THE ROUGH STATE 

We start our consideration from the high temperature 
region above the roughening transition when the crystal 
surface is in the atomicahy rough state. Since in this 
case the surface tension a has no singularity at the close- 
packed facets and weakly depends on the crystallographic 
direction, the tensor of surface stiffness is close to 
a^ifc. Thus, to simplify the analysis, we neglect any dis- 
tinction between surface tension and surface stiffness in 
the formulas for the Laplace pressure. 




FIG. 2: The isotropic growth of a "^He at 1.28 K. The diameter 
of the visible margin is 12 mm. 



A. Flat shape of the crystal surface 



B. Spherical shape of a solid 



Crystallization waves, representing small oscillations of 
the superfluid-solid ''He interface, are predicted by An- 
dreev and Parshin The spectrum of crystallization 
waves in the field of gravity can be found from the equa- 
tion [il 
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Here q is the wave vector, Ap = p' — p where p' and p 
are the densities of the solid and liquid phases, K is the 
interface growth coefficient, and g is the acceleration of 
gravity. The positive magnitude 5 > corresponds to 
the usual situation when the ^He crystal lies under the 
liquid phase. In this case the frequencies of small inter- 
face oscillations have the negative imaginary parts for all 
wave vectors and the crystal surface is always stable. 

Provided the solid phase occupies the half-space above 
the liquid, the dispersion equation remains the same but 
parameter g becomes negative. This means that the last 
term in ([T]), which dominates for sufficiently small g, can 
result in the positive imaginary part for the roots of the 
dispersion equation. The critical magnitude of the wave 
vector go is given by 
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where A '^l mm is the capillary length. The instability 
will develop faster for small wave vectors. For the har- 
monic q 0, the estimate for the time of developing the 
instability is given by 



Ap gK 



(2) 



For T ^1.1 K, the inverse growth coefhcient 1 /K ^ 2 m/s 
and t ^ 2 s. This is in a qualitative agreement with the 
observations of work Thus, the instability similar to 
the Rayleigh one develops at the lower surface of a ^He 
crystal in the atomically rough state. 



An interesting situation appears with the nucleation 
and growth of a crystal at the needle point at the tem- 
peratures above the roughening transitions (Fig. [2|) . In 
the hydrostatic pressure gradient the crystal grows in 
the shape of a ball which melts in the upper part and 
crystallizes in its lower part at the same time. In out- 
ward appearance this looks like the motion of a crystal 
downward. In work 0], in which the similar situation is 
studied, it is shown that the crystal with the isotropic 
growth coefficient conserves its spherical shape, and the 
descending motion velocity of the sphere v is equal to 



v^K^gR{t), 
P 



(3) 



where R{t) is a time-dependent radius of the sphere. 
Then, as follows from the center of the crystal shifts 
downward at a constant velocity in the hydrostatic pres- 
sure gradient provided the crystal volume remains un- 
changed. 

Let us consider stability of the spherical shape of a 
crystal against small shape perturbations in the field of 
gravity. For analysis, we choose the frame which origin 
is put at the center of a crystal. In this frame the cen- 
ter of a crystal is fixed and the liquid phase circulates 
around the crystal, outfiowing from the upper part of 
the sphere and flowing into the sphere in its lower part. 
We introduce ({t) ^ i? as a small perturbation of the 
radius R{t). The liquid phase is assumed to be incom- 
pressible and the solid one is motionless. In experiment 
(3| the velocities of the liquid flow are small and do not 
exceed 0.1 mm/s. For the hydrodynamic pressure pv'^/2, 
we have an estimate 10~^ dyne/cm^. This magnitude is 
much smaller than the typical hydrostatic pressure drop 
ApgR. In what follows, we neglect quadratic terms in 
velocity of the liquid phase. 

Next, in accordance with the experimental conditions, 
we take into account that the experimental cell is closed 
and no matter comes from outside. In other words, the 
total mass of the liquid and solid phases remains un- 
changed. In the frame comoving with the crystal the 
center of the crystal is fixed. As a result, we can omit 
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the spherical harmonic with I = from consideration 
since this harmonic is associated with the change of the 
volume. The next harmonic / = 1 is responsible for the 
displacement of the sphere as a whole and corresponds 
to the circulation of the liquid around of the sphere. We 
also omit this harmonic from our consideration. 

The determination of the oscillation spectrum is not 
difhcult but cumbersome. We give a scheme of solution 
and then the final result. Let axis z run in the vertical 
direction, and we seek for the general solution, expanding 
perturbation C in spherical harmonics. Since the liquid is 
assumed to be incompressible, it is convenient to describe 
its motion in terms of velocity potential (/> according to 
V = V</>. The mass flow across the interface is propor- 
tional to the chemical potential difference A/i 

J = p'KAfi. 

The continuity of the mass flow across the interface allows 
us to relate the velocity of the liquid at the interface with 
the interface growth rate (. After some calculations and 
involving the axial symmetry with respect to axis z, we 
arrive eventually at the dispersion relation 
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and ai is the perturbation amplitude corresponding to 
l-ih. spherical harmonic from a sum C,{t) — '^iO-i{t)Yi. 
Thus, we have a determinant of the inflnite order for de- 
termining the proper values. In the lack of gravity A = oo 
the off-diagonal terms vanish. The roots of the master 
equation give the oscillation spectrum of a spherical solid 
with damping. All the frequencies have a negative imag- 
inary part, resulting in the conclusion that in this case 
the spherically shaped solid is stable against small per- 
turbations of its equilibrium shape. 

In the presence of gravity we must involve the off- 
diagonal terms in Q. Note that the diagonal terms in- 
crease proportional to /'^ as / — )■ oo, while the off-diagonal 
terms remain finite and have the order of [R/X)^. This 
means that, as the degree of harmonic / increases, the 
relative effect of the pressure gradient reduces, agreeing 
with the stability of the flat interface at large wave vec- 
tors. As in the case of the flat interface, the instability 
develops in the first turn for the minimum degree of har- 
monics 1 — 2. 

The spherical shape of helium crystals in the hep phase 
is observed above the first roughening transition temper- 
ature at T >1.25 K (Fig. [5]). On the other hand, the hep 
phase is limited by the hcp-bcc transition at T =1.44 K. 




FIG. 3: The initial stage for the development of instability 
for a spherical crystal with the size larger than the critical 
radius. The dashed circle is the initial shape. 



Within this temperature range the kinetic growth coeffi- 
cient is small and varies insignificantly with the tempera- 
ture from 1/K ^11 m/s at the bcc transition to ~3 m/s 
near the roughening transition. Inserting these values 
into Eqs. (|31[S]) and solving numerically the master equa- 
tion for proper frequencies, we determine the critical ra- 
dius Rc =0.61 cm. The radius of crystals studied in Q 
does not exceed this magnitude. This may serve as an ex- 
planation that all the crystals grown in that experiment 
display stability of their spherical shape. 

The proper vectors in ^ are determined with accuracy 
to their sign. For R > Rc, this gives two possible cases 
for developing the instability shown in Fig. |3l In the 
first case at the bottom there appears a hogging in the 
upward direction. In the second case we see a formation 
of constriction. Note that the result refers only to the 
initial stage of developing instability. 



C. Crystals of the limited sizes 

Another factor which stabilizes the fiat surface is a 
finiteness of the crystal size. In this case the spectrum 
of crystallization waves becomes discrete. The instability 
will appear at the minimum wave vector. To estimate, we 
replace the hexagonal crystal shape with the equivalent 
cylindrical one of radius R. Then we can show that dis- 
persion equation ([l} holds its form but the perturbation 
amplitude C of the flat crystal surface is given by 



Cir, t) = CoJyn{qr)e 



where is the azimuth angle, J,„ is the Bessel function, 
and r is the distance from the center of a cylinder. For the 
circular symmetry at the fixed volume of a crystal, the 
minimum wave vector go is determined from the relation 



Ji((7oi?) = 0, (7o«3.83/i?. 



(6) 



From expressions (HHH) and ^ one can obtain that the 
instability should appear for a crystal with the diameter 
larger than 6 mm. 
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III. INSTABILITY OF A CRYSTAL FACET 



Two types of functions can be solutions of the equation 



Below the roughening transition temperature a singu- 
larity appears as a nonanalytic angular dependence in the 
behavior of the surface tension versus angle 9 between the 
direction of crystallographic axis and the normal to the 
crystal surface [1] 



a{e) = ao + P{T) I e I 
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Here (3 is the quantity proportional to the energy of a 
crystal facet step and is positive below the roughening 
transition temperature. This entails that the surface ten- 
sion and surface stiffness represent the drastically differ- 
ent quantities. The surface stiffness becomes infinite for 
zero angle 6 = Q. The singularity prevents from varying 
the shape of the crystal facet and qualitatively changes 
the conditions for stability. Let us consider the distinc- 
tions in appearing the Rayleigh instability in comparison 
with the case of the rough crystal surface. 



A. Rayleigh instability of the horizontal crystal 
facet 

Below the roughening transition the energy variation of 
the surface for sufficiently small perturbation amplitude 
C,{x, y) from the horizontal flat facet is equal to 



AS 



(VC)'+/3| vci 



1 ^ 



dx dy. 



The smallness of C and VC for the correct and consistent 
use of the expansion of energy in ^ is provided by in- 
equality /? <C OfQ. Let us start from the one-dimensional 
case, namely, flat bar of length L. The energy per unit 
length AE/D reads 
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IT 
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where D is the width of the crystal bar in the transverse 
direction. The boundary conditions at the ends of the 
bar 



C(-i/2) = C(V2) = 



(7) 



correspond to the case when the crystal surface is im- 
mobile at these points. It is convenient to introduce 
the dimensionless quantities according to a:' — x/L, 
rj = C{^pgL)/l3, and 7 = \/L. Then, 
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The extremum satisfles the equation 

j^ri" + 25{ri')f]" -I- 77 = 0. 
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'?(^')='?"|sin(x'-x[,)/7 



i.e., either a constant or a sine function. 

The height of a potential barrier for destructing the flat 
crystal facet as well as the type of a critical fluctuation 
depend on the magnitude 7. For sufficiently small length 
of the facet L < Ld = ttA, there exists a single trivial 
solution ri{x') — 0, and the flat crystal facet remains 
stable. For L > Ld, there appears a nontrivial solution 
which consists of two half-sinusoids and flat segment 
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Here the function ri{x') is continuous together with its 
derivative at points x' = ±x'q and vanishes at x' = 
1/2. The critical perturbation amplitude of fluctuation 
is equal to 



Vcl 
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or Cc 
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The height of the potential barrier, separating the tran- 
sition of the flat crystal facet to a distorted state, is given 
by 



AE 
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(8) 



As the size of a crystal facet increases, there appear ad- 
ditional possibilities for other fluctuations consisting of 
a combination of flat segments and half-sinusoids. The 
corresponding solutions, composed with n flat segments, 
appear as L > Lcn where Lcn = nLci- As the number of 
the flat segments increases, the perturbation amplitude 
of critical fluctuations does as well 
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The potential barrier height grows as the number n of 
possible flat segments increases 
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Provided the experimental conditions correspond to the 
conservation of the total mass including the both liquid 
and solid phases, the solutions should satisfy an addi- 
tional requirement 



L/2 



Cix) dx = 0. 

-L/2 

In this case the nontrivial solutions for critical fluctua- 
tions can be realized for the even numbers n and, corre- 
spondingly, flrst critical length becomes equal to Lc2- 
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For the temperatures well below the roughening tran- 
sition, the facet step energy coefficient /? is measured 
[H, and the numerical estimate of the coefficient in 
Eq. ([5]) gives the magnitude of the potential barrier about 
10^^ erg or 10^^ K for the basal facet of size ^1 cm. The 
overcoming of such barrier is practically impossible dur- 
ing the time of experiment. Thus, the appearance of the 
singular angle dependence in the function a{9) results 
in a drastic change for stability conditions of the crys- 
tal surface. With very small cooling below the roughen- 
ing transition temperature the potential barrier becomes 
sufficiently large and the probability of its overcoming 
due to thermal or quantum fluctuations is vanishingly 
small. This conclusion is confirmed by the experimen- 
tal evidence for stability of the crystal shape below the 
roughening transition. 

The one-dimensional problem and boundary condi- 
tions (O are chosen as simplest ones in order to illus- 
trate the method of solution and to obtain an analyti- 
cal estimate of the potential barrier height for develop- 
ing gravitaional Rayleigh-like instability. The full prob- 
lem should be solved employing the real crystal shape 
with the surfaces connecting the facets. We have ana- 
lyzed an onset of instability at the circular facet on the 
analogy with Sec. IIC. The mathematical treatment be- 
comes more complicated but the final result for the bar- 
rier height differs from Eq. (|5]) with a numerical coeffi- 
cient of about unity. 



B. Instability of the tilted crystal facet 

Here we analyze the stability of a crystal facet tilted 
with angle (jx jl — ip) against its small perturbations. The 
variation of energy Ai? per width of the crystal facet 
reads 



Ap g{xC,{x) GO'S, -\- 



sm ip' 
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where (^{x) is a perturbation amplitude taken from the 
flat surface. The angle (p = means the vertical position 
of the crystal and (p ~ Tr/2 corresponds to the horizontal 
position. 

Let axis Oz run along the normal to the crystal facet. 
The axis Oy is perpendicular to the axis Oz and the 
acceleration of gravity as well. The third axis Ox, lying 
in the plane of the facet, is perpendicular to the axes Oy 
and Oz. The axis Ox runs at the angle if to the direction 
of the acceleration of gravity. 

Let perturbation amplitude of the crystal facet C{x) 
from its initial position =0 take place in the direc- 
tion normal to the facet surface and be independent of the 
coordinate y directed horizontally along the crystal facet. 
In addition, we assume that the small perturbation am- 
plitude C{x) is finite within the region —L/2 < x < L/2 



and the following boundary conditions are fulfilled 

C(±i/2) = 0. (10) 

Here L is implied as a size of the crystal shape fiuctua- 
tion. The surface fiuctuations satisfy the conservation of 
the total mass of the liquid and solid phases 
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For ip = 7r/2, expression © goes over to the expression 
for a horizontal crystal facet with the small vertical per- 
turbations. For if = 0, we have the expression for the 
energy of a vertical crystal facet with the small pertur- 
bations in the direction normal to the facet. 

To find the minimum magnitude of the potential bar- 
rier which prevents from the development of instability, 
we must consider extremum of functional The ex- 
tremum of functional © satisfies the equation 

{aoC + l3 agnC')' + Apg{x cos ip + C sin ip) =0. (12) 

Provided the derivative C,'{^) does not vanish for all x 
within -L/2 < x < L/2, i.e., 

C{x) + 0, (13) 

equation (jl2p takes the simple form 

a^C^' {x) + Apg{xcosp) C,{x)smLp) =0. (14) 



If the condition (|13p is satisfied, the general solution of 
Eq. (HU reads 



C(a;) = —xcotp) -I- Asa\{x/\^p) -I- i?cos(a;/Ai^). (15) 

Unknown coefficients A and B are determined by condi- 
tions (fTO|) and ([Tl]) . Finally, we arrive at 
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Here A;^ plays role of an effective capillary length and 
determines the typical scale of length at which Q(x) varies 



A 



Y Ap g sin (p ^/sinp 

Note that, as — ^ for the vertical position of a crystal 
facet, the typical length diverges. 
Substituting (fTH]) into dH), we have 

AE ao A3 cosV f^/^ ( ^ 

= T- / rfa; a COS 2x — 2a COS X 

D 2 A^ swap J_i/2 V 

+ l + x^ + 2{P/aQ) I a cos a; - 1 I tan (/? ) (17) 



Here we have introduced the dimensionless quantities 
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The magnitude of energy p7|) can be written as 
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(18) 



which the facet step energy vanishes (Fig. |4]). In Fig. |4] 
the image of a ''He crystal at 0.901 K is given. Note 
that the lateral faceting disappears at 0.910 K, demon- 
strating the roughening transition for the a-facets. The 
crystal has a clear lateral facet which remained stable in 
the course of experiment during, at least, 10 min. 



For L = L/X^ ^ 1, equation (fTS)) gives 



IV. CONCLUSION 
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As a result, for sufficiently small-sized fluctuations with 
L <^ Lcr, the energy of fluctuations is positive AE > 
and such fluctuations are of low probability. The tilted 
facet is practically stable for such perturbations. The 
critical length is equal to 



, , ,5/3 4v/5/3/(3zApg) 
V oao y/cosLp 

This means that the crystal facet is practically stable if 
its size does not exceed the critical length. In fact, the 
magnitude of the energy barrier AE' ~ f3'^/Apg proves 
to be about 10~^ erg or 10^^ K. As the size of a crystal 
becomes larger than the critical length L > L^, the tilted 
crystal facet becomes unstable against distortion of its 
shape. 

The analysis of the crystal facet stability for the larger 
lengths L ^ X^p is simple. Let us represent Eq. ([1 
the form of the following inequality 



The gravitational instability at the atomically rough 
^He surface, which develops under the lack of any poten- 
tial barrier, is similar to the classical Rayleigh instability 




FIG. 4: The lateral facet of a *He crystal is stable. The tem- 
perature is below the second roughening transition by about 
20 mK. 



AE ^ aoL X^ cos^ ip 
13 ^ 4 sinp 



L 

L cot — 
2 



/3 ^ cosip 



L/2 



L/2 



dx (I a cos a; | +1) 



The right-hand side of the inequality is obviously nega- 
tive for L > TT, i.e. AE < 0. This means that a crystal 
facet distortion with such lengths is energetically favor- 
able. The tilted crystal facet becomes absolutely unstable 
if the facet size L > X^p. 

In conclusion, such high potential barriers can explain 
a gravitational stability of a crystal facet in the immedi- 
ate vicinity of the roughening transition temperature at 



when the heavier liquid lies above the lighter liquid. The 
distinction is that the time necessary for the develop- 
ment of instability is determined by the kinetic growth 
coefficient of a crystal surface. 

As for the smooth faceted crystal surface, having a 
singularity in the surface stiffness, the development of 
the surface instability becomes possible due to thermal 
or quantum fluctuations if the size of the facet surface 
exceeds the critical one. However, the large height of a 
potential barrier makes its overcoming impossible for a 
experimentally reasonable time. This explains the sta- 
bility of the lower facet for a free-growing ''He crystal 
observed in the experiment during a few hours. 
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